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Abstract 

A long X-ray flash was detected and localized by the instruments aboard the High Energy Transient 
Explorer II (HETE-2) at 00:03:30 UT on 2004 September 16. The position was reported to the GRB 
Coordinates Network (GCN) approximately 2 hours after the burst. This burst consists of two peaks 
separated by 200 s, with durations of 110 s and 60 s. We have analyzed the energy spectra of the 1st 
and 2nd peaks observed with the Wide Field X-Ray Monitor (WXM) and the French Gamma Telescope 
(FREGATE). We discuss the origin of the 2nd peak in terms of flux variabilities and timescales. We find 
that it is most likely part of the prompt emission, and is explained by the long-acting engine model. This 
feature is similar to some bright X-ray flares detected in the early afterglow phase of bursts observed by 
the Swift satellite. 
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1. Introduction 

X-ray flashes (XRFs) are generally thought to be a sub- 
class of gamma-ray bursts (GRBs). The main difference 
between XRFs and GRBs is the energy of the emission; 
the peak energy -E pea k of XRFs is distributed in the range 
from a few keV to 10—20 keV while the one of GRBs is 
distributed from -20 keV to -MeV (Barraud et al. 2003). 
Other properties such as timescales or features of the light 
curve are similar for XRFs and GRBs. Using the logarith- 
mic flucncc ratio log[5x/ <S 7 ] to categorize bursts, where 
S x is the 2-30 keV fluence and S 1 is the 30-400 keV 
flucncc, Sakamoto ct al. (2005) found that XRFs, X-ray 
rich GRBs (XRRs), and GRBs form a continuum in the 
[Sx, S^-plane and in the [Sx/S^, i? poa k]-plane. This is 
an evidence that all three kinds of bursts are the same 
phenomenon. Theoretical models which have been pro- 
posed to explain these soft events include off-axis view- 
ing (Yamazaki et al. 2002), a structured jet (Rossi et al. 
2002), and high-z GRBs (Heise et al. 2001). 

In this paper, we report the detection and localization 
of XRF 040916 by the HETE-2 satellite (Ricker et al. 
2002) and present the results of a detailed temporal and 
spectral analysis. As this burst has two peaks within a 
total time interval of —350 s, we discuss the origin of the 
long timescale. 

2. Observation 

2.1. Localization 

XRF 040916 triggered the WXM instrument on 2004 
September 16, at 00:03:30 UT (GPS: 779328222.72). This 
burst consists of two peaks lasting about 110 s and 60 
s separated by a time interval of —200 s. The HETE-2 
WXM instrument triggered on the 2nd peak. The ini- 
tial burst position was based on a rapid ground analysis 
using the WXM data, and was R.A. = 23 h 01 m 44 s , Dec. 
= -5°37'43" with a 90% confidence error circle of 18' 
radius. All coordinates in this paper are J2000. This 
position was reported in a GRB Coordinates Network 
(GCN) Position Notice at 02:26:16 UT (Yamazaki ct 
al. 2004) and is shown in Figure 1. A later, refined 
ground analysis using the WXM data gave an error box 
with the following corners: (R.A., Dec) = (23 h 02 m 01?68, 
-5°50'09"6), (23 h 00 m 23?76, -5°37'51f'6), (23 h 00 m 58f80, 
-5°19'51."6), (23 h 02 m 36?72, -5°32'09f'6) (90% confi- 
dence region). This was reported in a GCN Position 
Notice at 03:58:41 UT(Yamamoto et al. 2004), and is 
also shown in Figure 1. This error box is larger than 
those usually obtained by HETE-2 because XRF 040916 
was faint and the Soft X-ray Camera (SXC; 0.5—10 keV 
energy band; Villasenor et al. 2003) was not operating at 
the time. No other Interplanetary Network spacecraft ob- 
served this burst, so the localization could not be refined 
by triangulation. Identical localizations were obtained by 
using the data of the 1st and 2nd peak separately. 

The detection of the optical afterglow was first reported 
by Kosugi et al. (2004ab), who found it at R.A. = 
23 h 00 m 55?13, Dec. = -5°38'43."2 using SuprimeCam on 
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Fig. 1. HETE-2 WXM localization of XRF 040916. The 
circle is the initial 90% confidence region and the box is the 
refined 90% confidence region obtained by ground analysis. 
The point labeled "Subaru" is the location of the optical af- 
terglow (Kosugi et al. 2004, Henden 2004). 

the prime-focus of the Subaru 8.2m telescope (Figure 1). 
The afterglow was detected in the z', Ic, Rc, V and B- 
bands in this observation. Henden et al. (2004ab) also 
detected it in the Ic-band with the NOFS 1.0m telescope. 
The magnitude of the host galaxy was estimated to be 
fainter than Rc=25, which was the magnitude of the af- 
terglow measured 2 days after the burst. Despite these ob- 
servations, no redshift has been determined for this event. 

2.2. Temporal Properties 

Figure 2 shows the light curve of XRF 040916 in five 
WXM energy bands (2-5, 5-10, 10-17, 17-25, and 2-25 
keV). There are two peaks in the WXM bands and no 
significant emission above 17 keV. Table 1 gives the T50 
and T 90 durations in the 2-5, 5-10, 10-25, and 2-25 
keV energy bands for the first peak, the second peak, and 
the entire burst. 

The 1st peak is composed of two parts : a hard emission 
region and a soft one, referred to as regions (a) and (b), 
respectively, in figure 2. This 'hard-to-soft' evolution is 
typical of GRBs. The duration of the 2nd peak tends to 
be shorter at higher energies, which is a common feature 
observed in many GRBs ( Link et al. 1993; Fenimore et 
al. 1995). The duration of the 1st peak in the 10-17 keV 
band is longer than the one in 2—5 keV band because, in 
the high energy band the 1st peak consists of two pulses. 
In the region (a) there are few photons in the 2—5 keV 
band while there are certain photon contributions in the 
5-10 keV and 10-17 keV bands. 
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Fig. 2. Light curve of XRF 040916 in five WXM energy bands: 2-5, 5-10, 10-17, 17-25 keV (top to bottom). The light curve 
has been binned into 10 second bins. The regions (a) and (b) are the foregrounds of the 1st peak and (c) is the foreground of the 
2nd peak. Region (1) represents the background used for the 1st peak and region (2) is the background for the 2nd peak. 

Table 1. Temporal properties of XRF 040916 (Quoted errors correspond to la) 



Region 


1st & 2nd peak 


1st peak 


2nd peak 


Energy 
(keV) 


tgo tso 
(s) (s) 


tgo t.so 
(s) (s) 


tgo tso 
(s) (s) 


WXM 
2-5 
5-10 
10-17 


374.8 ± 6.2 287.5 ± 3.7 
350.2 ± 6.4 292.4 ± 4.1 
392.0 ± 17.3 325.6± 5.2 


103.3 ± 6.4 39.3 ± 3.2 
90.9 ± 1.7 50.4 ± 3.7 
120.5 ± 5.2 67.6 ± 4.6 


61.4 ± 4.1 28.3 ± 1.3 
39.4 ± 3.4 19.7 ± 2.1 
44.3 ± 8.6 18.4 ± 7.2 


FREGATE 
6-40 


349.0 ± 6.2 109.4 ± 5.4 


93.4 ± 13.1 43.0 ± 10.3 


55.3 ± 7.6 30.7 ± 9.4 



2.3. Emission between the 1st and 2nd peaks 

We calculated the count rates before and after the 1st 
peak, and compared them to the count rate after the 2nd 
peak. The results show that there is no significant differ- 
ence between them. We also tried to localize the burst 
using the photons in the time region between the 1st and 
2nd peaks. If there is emission from the burst between 
the peaks, the resulting localization should be the same 
as that of the 1st or 2nd peak alone. However, a significant 
localization could not be obtained, and we conclude that 
there is no significant emission between the peaks. The 
corresponding 2a upper limit is 1.8 x 10 -9 ergs cm -2 s _1 
(2-400 keV). 

2-4- Spectrum 

We performed spectral analyses of the two peaks sep- 
arately, and also of their sum. The background regions 
used were -460 < t < -330 s and -200 < t < -70 s for the 



1st peak, and -110 < t < -20 s and 90 < t < 220 s for the 
2nd peak, where t is the time since trigger. 

Two types of data sets (burst data and survey data) are 
provided by the HETE-2 WXM and FREGATE detectors 
(Atteia et al. 2003). The burst data are only available 
when the detector triggers on a burst, while the survey 
data are recorded whenever the HETE-2 satellite is oper- 
ating. The trigger occurred for the second peak of XRF 
040916, and consequently both sets of data were available 
for it, whereas only the survey data were obtained for the 
1st peak. 

The burst data include time-tagged photon data, while 
the survey data produce time-integrated (4.92 s) data for 
each wire in the WXM proportional counters. (The WXM 
instrument is composed of twelve 1-D position-sensitive 
proportional counters.) In order to improve the signal to 
noise ratio, and consequently the spectral analysis, we ap- 
plied a cut to the WXM photon time- and energy-tagged 
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data (TAG data), using only the photons from the pixels 
on the six wires of the X-detector and the two wires of the 
Y-detector that were illuminated by the burst. Moreover, 
as the gain is not uniform at the end of the wires (Shirasaki 
et al. 2003), we used only the photon counts that regis- 
tered in the center ± 50 mm region of the wires. These 
optimized TAG data were extracted and used for the 2nd 
peak of XRF 040916. 

We used the XSPEC v. 11.3.0 software package 
(Arnaud, K. 1996) to perform the spectral analysis. We 
simultaneously fit the WXM and FREGATE data with 
the four following functions; (1) blackbody function, (2) 
power-law function, (3) cutoff power-law function, and (4) 
Band function (Band et al. 1993). Table 2 shows the re- 
sults of the spectral analysis of each peak. 

Although the spectrum is soft, several facts argue 
against a type I X-ray burst (XRB) as the origin. First, 
the XRB emission mechanism is blackbody radiation, and 
we obtain a large x 2 except for the 1st peak (a) when we 
use a blackbody model to fit the data. Second, XRBs tend 
to be found in the Galactic plane or globular clusters, and 
they emit persistent X-rays; the Galactic latitude of XRF 
040916 is b = — 56°03', and there is no known persistent 
X-ray source or globular cluster at this latitude. Finally, 
the detection of an optical afterglow for this burst excludes 
an XRB interpretation. 

The analysis of the two peaks separately does not give 
the best fit. However, considering part (a) of the 1st peak, 
we obtain a statistically significant improvement using the 
cutoff power-law or the Band function, compared to the 
single power-law function, and the spectrum has a signif- 
icant break at -E pca k (Table 2). But in the other regions 
we cannot constrain the parameters by using the Band 
function. 

In the high energy band (above 10 keV), the photon 
index (3 tends to be less than -2, and smaller than the 
photon index a. E°^ k can be constrained to be below 10 
keV. But the band below 10 keV is near the low-energy 
threshold of WXM instrument (2 keV). If E°^ s ak is near the 
low-energy threshold, the spectrum will appear to follow a 
power-law, even if it is actually a Band function. We can 
constrain E°^ k using a constrained Band function. This 
can be done both for pure power-law spectra and power- 
law times exponential spectra with the required curvature 
in the detector energy range, but only the high-energy 
part of the Band function is allowed to produce a pure 
power-law spectrum. This is described in Sakamoto et al 
(2004). This method is applicable when the spectra of the 
burst have E°^l k near or below the low-energy threshold 
of the detector. 

Applying the constrained Band function model to each 
interval of XRF 040916, we obtained a constrained E°^l k . 
The results are shown in Figure 3 as the posterior prob- 
ability density distribution for E°^ k . From these distri- 
butions, we find best-fit values for 68%, 95% and 99.7% 
probabilities (Table 4). We conclude that XRF 040916 
is extremely soft compared to typical GRBs. E P e a k was 
determined for the 1st peak, but not for the 2nd (it is less 
than 4.8 keV with 99.7% probability); this hard-to-soft 



evolution is common in GRBs. 

We also calculated the hardness ratio of the 30—400 keV 
(S-y) and 2—30 keV (S'x)flucnccs. Using this parameter, 
we can categorize bursts as XRFs when logfSx/fi^] > 0, X- 
ray rich GRBs when —0.5 < log[5x/<S 7 ] < 0, and classical 
hard spectrum GRBs when loglSx/S^} < —0.5. For the 1st 
peak, as log[5x/<S 7 ] < with the cutoff power-law func- 
tion, this burst would be classified as XRR rather than an 
XRF. For the 2nd peak, log[5 x /5 7 ] is 0.28, with a 90% 
confidence lower limit of 0.13, that is, the emission of the 
2nd peak is softer than that of the 1st peak. Considering 
finally the total emission of the 1st and the 2nd peaks, 
the value of log[5x/5 7 ], using the cutoff power-law func- 
tion, is 0.17 with a 90% confidence lower limit of 0.06. 
Therefore we regard the entire burst as an XRF. 

We consider the 1st peak to be clearly the prompt emis- 
sion of XRF 040916; we discuss the nature of the 2nd peak 
in the next section. 

3. Discussion 

3.1. Redshift estimates 

The spectral softness could be explained if XRF 040916 
were a high-redshift GRB. To check this hypothesis, 
we first estimated the redshift from the Amati relation 
(Amati et al. 2002) using the 1st peak. Amati et al. 
(2002) found a relation between the isotropic-equivalent 
radiated energy E lso and the burst-averaged value of Speak 
in the rest frame (E^ k tx E? s £, with E r p c ^ k in keV and 
E iso in units of 10 52 ergs). 

Assuming this relation, it is possible to estimate the 
redshift using only the flux and S pea k of a burst. We 
calculated E^£ k and Ei so from the spectral parameters 
of the 1st peak assuming various redshifts. As shown in 
Figure 4, the smaller the redshift is, the more consistent 
the computed values are with the Amati relation. This re- 
sult is also consistent with the Subaru redshift constraint 
of z < 3, imposed by its detection of the optical afterglow 
in the B-band(Kosugi et al. 2004) 

We have also computed an upper limit to the pseudo- 
rcdshift (with the method described in Pelangeon et al. 
2006) using the WXM spectrum for the most intense 15s 
long part of the 1st peak. To do this, we derived the upper 
limit to Epeak using the constrained Band model to fit the 
data. We find that Speak < 6.2 keV with 90% confidence, 
leading to a pseudo-redshift < 0.7, consistent with the 
preceding results. Thus we can reject a high redshift for 
this burst. 

3.2. Origin of the 2nd peak 

In the spectral analysis, we could not determine Speak 
for the 2nd peak, but we found that above 10 keV the 
best fitting model was a power-law, and that S pea k lay 
in the low energy region (<4keV). Furthermore, the 1st 
and 2nd peaks are separated by ~ 200 s. This suggests 
that the 2nd peak could be the X-ray afterglow rather 
than prompt emission. Indeed there are at least five 
possibilities to explain it. (1) beginning of the afterglow 
(2) Ambient density fluctuations of the X-ray afterglow, 
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Table 2. Results of the spectral analyses performed for XRF 040916 



(Note. The quoted errors correspond to 90% confidence. ) 



Region 



Function 



kT 
[keV] 



rpo bs 
peak 

[keV] 



X 

(d.o.f.) 



1st peak (a) blackbody 5.4 

Power law 
Power law (above 10 keV) 
Cutoff power law 
Band 

1st peak (b) blackbody 1.2 

Power law 
Power law (above 10 keV) 
Cutoff power law 

1st peak (a+b) blackbody 

Power law 
Power law (above 10 keV) 
Cutoff power law 

2nd peak blackbody 2.4 

Power law 
Power law (above 10 keV) 
Cutoff power law 

1st + 2nd blackbody 

Power law 
Power law (above 10 keV) 
Cutoff power law 



TTJT 

-1.2 



+0.3 
-0.2 



o q+0.8 
z -°-0.6 



+0.9 



+0.4 

-0.4 



-1 fi +0 - 2 

-2 +0 - 3 

Z ' U -0.5 
}+0.39 
-0.95 

-0 lfi+ 2 08 
U - 1D -1.03 



-0.39: 



9 9+0.4 
^•^-0.3 
9 q+0.7 
— 1.7 
_ X g+0.7 a 



2+0.1 
-0.2 
vK>.3 
-0.5 
i 7+0.3 
_i -'-0.2 



-li 

-2.0 H 



^+0■2 
-0.3 
;+0.6 
-1.3 

-1 6+ 08 

± - D — 0.5 



-1.9 H 
-2.5 H 



.1 q+o-i 

9 9+O.3 
z ' z -0.5 
1 S+0-3 
j -0.3 



-2 4 +0 - 7 



90 9+48.2 

25.6t?°j 



+3.7 a 



3.8 



70 q+°° 
' a -°-59.5 



99 q+CJO 
zz -°-18.3 



25.8l:^. 5 



21.68(21) 

24.41(21) 

8.76(13) 

19.19(20) 

18.00(19) 

25.57(21) 

18.02(21) 

13.37(14) 

18.37(20) 

55.60(34) 

26.00(34) 

16.69(22) 

25.66(33) 

37.80(17) 

19.53(17) 

10.13(10) 

18.86(16) 

82.30(47) 

42.16(47) 

19.45(28) 

41.17(46) 



The data do not allow a determination of the lower limits for the cutoff power-law function a or E°^l k parameters 



Table 3. Flucnccs using the cutoff power-law function 



region 


Sx 


S-y 


log[5 x /5 7 ] (90% lower limit) 




[10~ 7 erg cm~ 2 ] 


[10 -7 erg cm -2 ] 




1st peak 


4.92±0.72 


6.18±0.90 


-0.10 (>-0.23) 


2nd peak 


2.91±0.50 


1.53±0.26 


0.28 (>0.13) 


1st peak + 2nd peak 


7.74±0.91 


5.27±0.62 


0.17 (>0.06) 



Table 4. Results of posterior probability density 





best-fit value (keV) 


£ p °e b ik (keV) 


£ p °c b ik (keV) 


^e b Ik (keV) 






with 68% prob. 


with 95% prob. 


with 99.7% prob. 


1st peak (b) 


3.8 


1.0 - 4.9 


<5.9 


<7.1 


1st peak (a) + (b) 


3.9 


2.5 - 4.4 


<4.8 


<5.5 


2nd peak 




<2.7 


<3.6 


<4.2 


1st peak + 2nd peak 




<2.3 


<3.0 


<3.5 




Fig. 3. Posterior probability density distribution for E^^,. The vertical solid lines define the 68% probatility interval for E** k , 
while the dashed and dotted lines show the 95% and 99.7% probability upper limits on E°^ k . Figures from top left to bottom right 
show the regions for the 1st peak (b), the 1st peak (a) + (b), the 2nd peak, and the 1st + 2nd peaks, respectively. 
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BeppoSAX-GKBs (Amati et al. 2002) and HETE-2 GRBs. The dashed line is the equation, E T p e ^ = 89(£ iso /10 52 ergs) ' 5 kcV 
given by Amati et al. (2002). 



No. ] 

(3) Patchy shell, (4) Refreshed shock and (5) Long-acting 
engine model (loka et al. 2005). We consider them in the 
following sections. 

3.2.1. Beginning of the afterglow 

Piro et al. (2005) said that GRB 011121 and GRB 
011211 showed a late X-ray burst taking place a few hun- 
dreds of seconds after the prompt emission. The spectral 
and temporal evolution of the afterglow indicated that 
fireball evolution took place in the interstellar medium 
for GRB 011211 and in the wind for GRB 011121 respec- 
tively. In both cases, the decay time of a late X-ray burst 
(i.e. afterglow) is supposed to be proportional to ~ t^ 1 
(Sari et al. 1998). 

If we set the zero time (to) of the emission episode 
for the 1st peak, then the decay time of the 2nd peak 
is proportional to t~~ 12J±0 - 3 in 2—25 keV energy band. 
Since this result is inconsistent with the above afterglow 
model, we can reject the possibility of the beginning of 
the afterglow. 
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3.2.3. Patchy shell model 

In the patchy shell model the variability timcscale of 
the afterglow at time t must be At > t (Nakar & Oren 
2004). In this model, the GRB jet consists of many 
subjets (Yamazaki et al. 2004; loka & Nakamura 2001). 
Since we observe an angular size 7 -1 within a GRB jet 
with Lorentz factor 7, the flux depends on the angular 
structure for the observer. The patchy shell model cannot 
make a bump with variability timescale At ~ 60 s, which 
is shorter than the observed timescale t ~ 370 s, so we 
can marginally reject this model from the point of view 
of the variability timescale. In addition, the flux of the 
2nd peak is as bright as that of the 1st one. If this burst 
is explained by the patchy shell model, we must assume 
very large shell non-uniformity. In this case, we have to 
assume that we first observe the dark part of the shell, 
and then the bright part of the shell as the Lorentz factor 
drops. As the assumption of a large non-uniformity shell 
is unrealistic, we can consequently exclude the patchy 
model. 
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3.2.2. Ambient density fluctuation model 

This model explains afterglow variabilities by ambient 
density fluctuations caused by turbulence in the interstel- 
lar medium or variable winds from the progenitor star. If 
XRF 040916 was in the afterglow stage when the 2nd peak 
occurred, we can use the following formula for a kincmat- 
ical upper limit on the variability (loka et al. 2005) 



F u ~ Jc vF v 



t) 



(1) 



where At is the variability timescale, t is the observed 
time since the burst, F v is the assumed power- law 
baseline of the afterglow flux, F is the bolometric base 
flux, f c ~ (v m j V c )( p-2 )/ 2 where the cooling frequency is 
v c , the characteristic synchrotron frequency is v m , and 
the electron power-law distribution index is p, F jvF v ~ 
(y I V c )( p ~ 2 )/ 2 , and At, AF V are the timescale and am- 
plitude deviations above the baseline, respectively. This 
formula gives the maximum possible afterglow variability 
due to ambient density fluctuations. We can estimate 
the factor F/vF^ assuming the standard afterglow model 
for v m <v c <v (the X-ray band at t ~ 370s). Since v c 
~ 10 16 Hz at t - 370s for p ~ 2.2 (Sari et al. 1998), 
we have F/vF v ~ 1 for the X-ray band (y ~ 10 19 Hz). 
Furthermore, substituting At ~ 60 s in equation 1, the 
right-hand side becomes ~ 0.1. But we can clearly see 
that AFv/F u ;§> 1 from the flux variabilities in Figure 
2, because the flux Fv in the region between the 1st 
and 2nd peaks (< 1.8 x 10 -9 ergs cm~ 2 s" 1 with 2a 
confidence level, 2—400 kcV) which is assumed to be the 
afterglow stage, is less than the flux of the 2nd peak 
(8.3 x 10" 9 ergs cm" 2 s - \ 2-400 keV) by a factor of 
more than 5 with 2a confidence level. Since the observed 
variability exceeds the maximum allowed value, ambient 
density fluctuations cannot explain the nature of the 2nd 
peak. 



3.2.4- Refreshed shock model 

In this model, multiple shells are ejected at various ve- 
locities, and the variability occurs when the slow inner 
shell catches up with the fast outer shell a long time later, 
since the velocity of the outer shell decreases through the 
interaction with the ambient medium (Rees & Meszaros 
1998; Panaitescu et al. 1998; Kumar & Piran 2000; Sari 
& Meszaros 2000; Zhang & Meszaros 2002). The variabil- 
ity timescale at time t is given by A t > t/4 (loka et al. 
2005). In the case of XRF 040916, however, the variabil- 
ity timescale is A t < t/A for t ~ 370 s and At ~ 60 s, so 
we can marginally reject the refreshed shock model. 

In addition, the flux variability will not be equal to 
the GRB flux. The relation between the energy increase 
factor F and the flux increase / is given by / = F^ 3+p ^ 4 
for v m < v < u c , and / = F( 2+p )/ 4 for u m ,p c < v (Granot 
et al. 2003). In either case, / ~ F with typical values of 
p ~ 2.2. Because / 3> 1 from the observation, we must 
have F 3> 1. This means that the slow shells have very 
large energy, which contradicts the refreshed shock model 
(Refreshed shocks produce changes of less than 1 order 
of magnitude). Therefore this model, too, can be rejected. 

3.2.5. Long- acting engine model 

In this model, at the observed time t, the central en- 
gine is still active and emitting shells (Rees & Meszaros 
2000; Zhang & Meszaros 2002; Dai & Lu 1998). This can 
explain variability timescales down to a millisecond and 
there is no restriction on the flux variability. The most 
likely explanation of the 2nd peak is therefore the long- 
acting engine model. Only this model can explain both the 
variability timescale and the flux variability for the 2nd 
peak. In this scenario, both the 1st and 2nd peaks have 
the same mechanism, and both must show the same spec- 
tral and temporal features as GRBs. As we have shown 
in Section 2.2, the timescale of the 2nd peak in Table 1 
is shorter at higher energies, which is a typical feature of 
GRBs. 
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Furthermore, we considered the curvature effect 
(Kumar & Panaitescu 2000; Liang et al. 2006). This 
effect is the rapid decay due to the observed receiving the 
progressively delayed emission from higher latitudes. The 
following formula represents the curvature effect, 

\-(l+/3) 

JM,) - it (^r) • (2) 

where /3 is the X-ray photon index during the decay, to 
is the time zero point of the emission episode related to 
decay and A is normalization parameter for decay com- 
ponent. We applied this effect to the 2nd peak. Because 
of the poor statistics in the 2nd peak, we used the X-ray 
photon index (/3 ~ 1 .9) during the entire 2nd peak (rise 
and decay) with the background subtracted in the region 
(2). From fitting analysis, we obtained to = -34.8 ± 34.9 
s, near the time zero point for the 2nd peak. This result 
implies that the zero time for the 1st peak doesn't coin- 
cide with the 2nd peak one: then the central engine was 
reactivated. 

And as we showed in Section 2.4, the spectrum of the 
2nd peak is probably softer than that of the 1st. Thus 
-Epeak appears to decrease with time and flux, in accord 
with the well known hard-to-soft evolution in GRB spectra 
(Fishman & Meegan 1995). This feature is also consistent 
with the long-acting engine model, and we consider it to 
be the most reasonable mechanism to explain the soft 2nd 
peak of XRF 040916. 

Furthermore, we can compare this event with the Swift 
bursts XRF 050406 and GRB 050502B. Swift (Gchrcls ct 
al. 2004) can detect early X-ray afterglows with the X- 
Ray Telescope (XRT; Burrows et al. 2005a) and it has 
detected bright X-ray flares 100—1000 seconds after the 
prompt emission for these bursts. They have been ex- 
plained by the long-acting engine model (Burrows et al. 
2005b, Romano et al. 2006, Falcone ct al. 2006). Because 
the XRT is only sensitive to photons in the 0.2—10 keV 
energy range, the spectral features of the X-ray flares are 
not as well resolved as those of the prompt emission, and 
their emission mechanism is still under discussion. In 
GRB 011211( Jakobsson et al. 2004; Holland et al. 2002; 
Piro et al. 2005), GRB 011121 (Piro et al. 2005) and GRB 
021004( Bersier et al. 2003; Halpern et al. 2002) the after- 
glows included X-ray flares, whose study revealed a wealth 
of information about the central engine and its surround- 
ing regions. For XRF 050406, the bright flare occurred 
^200 s after the prompt emission, and for GRB 050502B, 
it occurred ~700 s later . The hard to soft count rate ra- 
tios for the flares was similar to the spectral evolution of 
the prompt GRB emission. As the time interval of XRF 
040916 is also of the order of ^200 s and the emission of 
the 1st peak is harder than that of the 2nd according to 
the \og[Sx/ S-y) ratio, the timescale and the spectral evolu- 
tion between the peaks for XRF 040916 is similar to those 
of XRF 050406 and GRB 050502B, and can be explained 
by the long-acting engine model. 



4. Conclusion 

In this paper we have reported the HETE-2 
WXM/FREGATE observations of the XRF 040916, 
which consists of two peaks separated by ~200 s. For 
this burst, only an upper limit to the redshift has been 
reported based on optical observations (z < 3). Taking 
into account this limit, our estimate from the Amati 
relation is consistent with a small redshift. In terms 
of the spectral evolution, it seems that the 2nd peak is 
softer than the 1st peak. We have studied the emission 
mechanism of the 2nd peak and discussed its most 
probable origin. Considering different models ( beginning 
of the afterglow, ambient density fluctuations, patchy 
shell, and refreshed shock model), we found that the 
long-acting engine model is the most plausible one to 
explain both the timescale and flux variabilities. In some 
GRBs or XRFs (e.g. GRB050502B and XRF 050406) 
(Burrows ct al. 2005b), bright X-ray flares were observed 
and these bursts are also explained by the long-acting 
engine model. Consequently, our results indicate that the 
case of XRF 040916 is similar to those of the X-ray flares 
detected by Swift. 
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